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Section  1 


INTRODUCTION 


In  fluid-dynamics  experiments  the  need  to  provide  a  smooth  transition 
from  a  round  cross  section  to  a  square  one  is  often  encountered.  Recent 
advances  in  experimental  and  computational  studies  have  shown  that  the  shape  of 
the  transition  section  dramatically  impacts  the  downstream  flow  conditions.  The 
Advanced  Propulsion  Division  of  the  Aero  Propulsion  and  Power  Directorate 
(WL/POPT),  Wright-Patterson  Air  Force  Base,  funded  a  task  to  design  such  a 
transition  section  for  a  water  tunnel.  M  The  resulting  section,  although 
successful,  was  for  a  very  specific  design  parameter. 

The  purpose  of  this  report  is  to  provide  1)  a  generalized  formula  which 
will  generate  a  transition  surface  for  any  given  diameter  of  the  round  end,  the 
width  of  the  square  end,  and  the  length  of  the  transition  section  and  2)  a 
computer  code  which  will  generate  the  surface  coordinates  in  a  format  suitable 
for  direct  input  to  the  front-end  VAX  computer  currently  employed  in  the  Air 
Force  Machine  Shop  in  Bldg.  5  at  Wright-Patterson  Air  Force  Base.  This 
automatic  procedure  benefits  the  design  task  of  the  test  sections,  the  inlet 
diffusers,  and  the  converging  part  of  the  converging-diverging  section  of  a 
supersonic  nozzle. 


Section  2 


ANALYTICAL  FORMULATIONS 


The  criteria  for  designing  the  transition  section  are:  1)  a  zero  slope  in  the 
longitudinal  direction  at  either  end  of  the  nozzle  and  2)  an  inflection  point  along 
the  length  of  the  transition  section  at  a  specified  location.  Generally,  the 
maximum  inflection  occurs  about  two-thirds  of  the  way  down  the  transition 
section. 

& 

AXIAL  CUT  * 


» 


The  coordinate  system  chosen  for  the  analysis  is  displayed  in  Fig.  1 .  The 
axial  direction  is  denoted  by  x,  the  distance  from  a  surface  point  to  the  x-axis  is 
denoted  by  R,  which  is  a  function  of  x  and  the  tangential  angle  (0)  between  the 
point  and  the  x-z  plane.  The  width  of  the  square  end  is  denoted  by  SW1,  the 
diameter  of  the  round  end  by  RD1,  and  the  length  of  the  transition  section  by 
XLN.  Based  upon  these  notations  and  criteria,  a  formula  for  R  as  a  function  of 
x  and  0  was  given  [1]  in  non-dimensional  form  as  follows: 

F  =  80.4375  +  x*i]  (1 ) 

l4.5  2.75  6.5 

where  F  *  R  /  SW1  /  VI,  x  =  x  /  XLN. 
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Figure  1.  Coordinate  System  Used  in  Analysis, 
x  -  Downstream  Direction. 


CROSS-SECTIONAL  CUT 


The  shape  of  the  cross  section  of  the  design  transition  will  vary  from 
round  to  square.  The  general  form  describing  this  cut  is: 


y  jn(x) 


<*) 


fe'® 


1. 


(2) 


Since  B(x)  must  fulfill  Eq.  (1)  and  two  end-point  conditions,  it  must  follow  the 
form: 


B(x)  =  BK1  [xii-i£i  +  x^L]  +  BK2 
l4.5  2.75  6.5 

where  BK1  =  40.2187  *  (RD1  -  SW1), 
and  BK2  =  SW1  /  2.0 
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In  obtaining  the  expression  for  n(x),  another  condition  must  be 


1 


considered— the  axial  line  going  through  the  comer  of  the  square  section.  Since 
this  line  is  in  a  symmetric  plane,  the  following  condition  is  satisfied: 


i 


y  -  zs  (4) 

Again,  this  line  must  fulfill  Eq.  (1)  and  two  end-point  conditions;  thus, 

Z«(x)  =  ZK1  [*--5  - *-5-  +  +  ZK2  ,,, 

s  4.5  2.75  6.5  (5) 

where  ZK1  =  80.4375  (BJD1 .  SWL)  /  fY. 

2.  (2. 


and  ZK2  =  SW1 

2. 

> 


Combining  Eqs.  (5)  and  (2)  yields: 


n(x)  = 


•og(^) 

log(^) 

BOc)' 


(6) 


Since  B(x)  and  n(x)  in  Eq.  (2)  have  been  derived,  the  following  relation 
can  be  utilized  to  calculate  the  surface:  10 


R(x,9)  = 


B(x)(  l. 
cosO  1 .  +  tann^0 


(7) 


» 
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Then,  the  physical  coordinates  y  and  z  can  be  obtained  as: 


y  =  R  sin(0),  z  =  R  cos(0)  (8) 

Determination  of  the  four  constants  (BK1,  BK2,  ZK1,  ZK2)  in  Eqs.  (3) 
and  (5)  is  dependent  upon  ihe  geometric  constraints.  In  this  case,  the  round  end 
is  larger  than  the  square  one.  For  the  reverse  condition,  the  values  of  those 
constants  must  be  reformulated  as  follows. 

BK1  =  80.4375  *  (SW1  -  RD1)  /  2. 

BK2  =  RD1  /  2 

ZK1  =  80.4375  *  (SW1  /  2.  -  RD1  /  2.  /  V2) 

ZK.2  =  RD1  12.1  ii 
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COMPUTER  CODE  AND  DATA  FORMAT 


A  computer  code  which  yields  the  surface  coordinate  of  the  transition 
section  was  developed  and  is  listed  in  the  Appendix.  The  input  parameters  are 
as  follows: 


RD1 :  length  of  the  diameter  of  the  round  end 

SW1 :  width  of  the  square  end 

XLN:  length  of  the  transition  section 

NS:  number  of  axial  cuts 

NC:  number  of  cross-sectional  cuts 


The  units  for  the  diameter,  the  width,  and  the  length  should  be  the  same, 
e.g.,  inches.  The  choice  of  NS  and  NC  will  depend  upon  the  precision  required. 
Assigning  higher  values  to  NS  and  NC  will  result  in  a  finer  grid.  For  example, 
there  is  a  need  to  design  a  6-inch-long  transition  section  to  interface  a  rounded 
settling  chamber  having  a  diameter  of  5.875  inch  and  a  square  test  section 
having  a  width  of  1  inch.  Values  of  31  and  15  were  assigned  to  NS  and  NC, 
respectively.  The  computation  process  takes  advantage  of  geometric  symmetry. 
Only  one-eighth  of  the  entire  circumference  was  actually  computed.  Due  to  the 
symmetry,  the  output  yields  the  grid  of  an  entire  quarter.  Thus,  in  this  quarter, 
the  number  of  the  axial  cuts  is  two  times  NC,  i.e.,  30. 
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The  output  of  the  computer  code  is  the  grid  positions  which  are  the 
interception  points  of  the  axial  and  cross-sectional  cuts.  An  automatic 
machining  device  used  in  the  Air  Force  Machine  Shop  reads  in  the  grid  point  of 
each  individual  cross  section  following  the  ASCII  format: 

xl,  yl,  zl 
x2,  y2,  z2 


xn,  yn,  zn  A 

At  the  end  of  a  cross  section,  the  last  grid  coordinate  is  followed  by  a  space  and 
a  capital  letter  A.  This  format  is  used  repeatedly  for  the  remainder  of  the 
cross-sectional  cuts. 

The  three-dimensional  plot  of  the  grid  generated  for  a  quarter  of  the 
transition  section  described  in  the  above  example  is  displayed  in  Fig.  2.  The 
grid  coordinates  were  stored  on  magnetic  tape  at  a  density  of  1600  BPI  using  a 
non-compressed  format  and  given  to  the  shop.  The  completed  transition  section 
were  installed  in  a  supersonic  combustion  testing  facility. 


0.0 
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SUMMARY 


A  computer  code  was  developed  for  the  design  of  a  transition  section 
interfacing  round  and  square  sections  of  a  wind  tunnel— a  design  challenge 
which  is  encountered  constantly.  The  output  data  have  been  formatted  to  be 
readable  to  the  computer  currently  in  use  in  the  Air  Force  Machine  Shop  at 
Wright-Patterson  Air  Force  Base.  The  procedure  involved  in  converting  design 
to  product  has  been  successfully  tested  in  association  with  the  personnel  at  the 
Air  Force  Machine  Shop.  Since  the  program  was  written  in  standard 
FORTRAN  language,  it  can  be  run  on  any  computer.  The  output  format  can  be 
changed  according  to  the  requirements  of  different  facilities  and  presented  using 
the  available  graphics  package. 
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APPENDIX 


COMPUTER  PROGRAM  FOR  GENERATING 
SURFACE  OF  TRANSITION  SECTION 


PROGRAM  SRNZ3N 

(>»  SQUARE-ROUND-NOZZLETRANSmON 

C>» 

C  Drs.  Tzong  H.  Chen  &  Abdollah  S.  Nejad 

C  Aero  Propulsion  &  Power  Directorate 

C  Wright  Laboratory 

C  Wright-Patterson  Air  Force  Base,  Ohio 

C 

c - 

REAL*8  RX  ATI  T1 

DIMENSION  YF(300),  XL(300),  BX(300),  ANX(300), 

#  YX(300),  RX(300),  ZS(300),  AX(300),  ZC(300),  ZF(300) 

DIMENSION  X(100),  Y(25),  Z(100,25),  ZW(  1 00),ZR(300),  YR(300) 
DIMENSION  ICA(300) 

COMMON  /CNST1/  RD1,  SW1 

OPEN(3,FILE-*NOZZLE.XYZ',  STATUS  =  'UNKNOWN*) 

333  CONTINUE 
WRITE(*,701) 

701  F0RMAT(1X,'INPUT  DIAM.  (RND.  END),  WDTH  (SQUARE  END)  ■) 
READ(*  •)  RD1,  SW1 

WRITE(*,702) 

702  F0RMAT(1X,'INPUT  TUNNEL  LENGTH  *) 

READ(*,*)  XLN 

XMIN  =  0.0 
YMIN  =  0.0 
XMAX  =  RDl/2. 

YMAX  =  RDl/2. 

WRITER  *,355)  XMAX,  YMAX 
355  FORM AT(  1  X,'XM AX,  YMAX  ’,2(F10.4,2X)) 

C - 

BK1  =  80.43750001  *  (RD1-SW1)  /  2. 

BK2  =  SW1  /2.0 

C - 

ZK2  =  SW1  / 2  0 

ZK1  =  80.43750001  *  (RD1/2.  -  SW1/2.  * 

#SQRT(2.0))/SQRT(2.0) 

WRITE(*,703) 

703  FORMAT(  1  X,'BK  1 ,  BK2,  ZK 1 ,  ZK2  ') 

WRITE(*,950)  BK1,BK2^K1,ZK2 

950  FORM  AT(  1  X,4(F  1 0.4, 1 X)) 

WRITE(*,704) 

704  FORMAT(  1 X, ’INPUT  NS:  NO.  OF  CROSS  SECTION  ’) 

READ(V)NS 

C  NS  ®  25 
C  NC  =  31 
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WRITE(*,705) 

FORMAT(  IX, 'INPUT  NO.  OF  AXIAL  CUT 
READ(V)  NC 
06  =  1.  /(NS  -  1) 

CL1  =  ALOG10(0.5) 

DO  60  I  =  1,  NS 

XL(I)  =  C16  *  FLOAT(I-l) 

BX(I)  =  CBX(BK1,BK2,XL(I)) 

ZS(I)  =  CZS(ZK1,ZK2,XL(I)) 


BB  =  BX(I) 

IF(BX(I).EQ.O.)  BB=  l.E-20 
AAL  =  ZS(I)  /  BX(I) 
IF(AAL.GT.O.)  GO  TO  616 
AAL  =  ABS(AAL) 


616  ANX(I)  =  CL1  /  ALOGIO(AAL) 

WRITER, 955)  XL(I),  ANX(I),  BX(I),  ZS(I) 
955  FORM  AT(  1  X,4(F  1 0.4, 1 X)) 

60  CONTINUE 


NC2  =  2*  NC  -  1 

DAI  =3.14159265  /  4./  (NC  -  1) 

DO  70  I  =  1,  NS 
DO  80  K.  =  1,  NC 
AX(K)  =  (K  -  1)  *  DAI 
Cl  =  COS(AX(K))  +  l.E-22 
T1  =  TAN(AX(K))  **  ANX(I) 
IF(T1.GT.1.E22)T1  =  1.E22 
ATI  =  1.  /  (1.  +  Tl) 

RX(K)  =  BX(I)  /  Cl  *  (ATI)  **  (l./ANX(I)) 
YX(K)  =  RX(K)  *  S1N(AX(K)) 

ZC(K)  =  RX(K)  *  COS(AX(K)) 

CONTINUE 

FORM  AT(  1  X,2(F8 .3,1  X)) 

CONTINUE 
IL  =  I 

DO  133  IN  =  1,  NC 
YF(IN)  =  YX(IN) 

ZF(IN)  =  ZC(IN) 

CONTINUE 
NCP1  =  NC  +  1 
DO  155  IN  =  NCP1,  NC2 
ID  =  IN  -  NC 
IA  =  NC  -  ID 
YF(IN)  =  ZC(IA) 

ZF(IN)  =  YX(IA) 

CONTINUE 
NP  =  NC2 

XLL  =  XL(IL)  *  XLN 
DO  273  HA  =  1,NC2 
ICA(IIA)  = '  ' 

CONTINUE 
ICA(NC2)  = '  A' 


i  < 
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0»  REVERSE  ORDER  TO  D  0  INTERPOLATION 
DO  177  IQ  =  1,NC2 
IV  =  NC2  -  IQ  +  1 
ZR(IV)  =  ZF(IQ) 

YR(IV)  =  YF(IQ) 

WRITE(3,535)  XLL,  ZF(IQ),  YF(IQ),  ICA(IQ) 
WRITE(*,535)  XLL,  ZF(IQ),  YF(IQ),  ICA(IQ) 

535  FORMAT(  1  X,F  1 0.4, 1  X,F  1 0.4, 1  X,F  1 0.4.A2) 

177  CONTINUE 
ENDFILE  3 
70  CONTINUE 
WRITE(*,735) 

735  FORMAT( IX,' MORE  COMPUTATION  I-YES.O-NO  ’) 
READ(*,740)  IMORE 
740  FORMAT(Al) 

IF(IMORE.EQ.  1)  GO  TO  333 

STOP 

END 

£**********************•**************♦***** 

FUNCTION  CF(X) 

AN  =  4.5 
BN  =  2.75 
CN  =  6.5 
AP  =  4.5 
BP  =  5.5 
CP  =  6.5 
A  =  X**AP  /  AN 
B  =  X**BP  /  BN 
C  =  X**CP  /  CN 
CF  =  (A-B  +  C) 

RETURN 

END 

c***„****,**,,,**,****,,*,***«*,,,*,********»»»* 

FUNCTION  CZS(CN1,CN2,X) 

C»>  EQ.  6 

CZS  =  CN1  *  CF(X)  +  CN2 

RETURN 

END 

C**M*MMmOM*»MOm**«*M»t***M*M»M*M 

FUNCTION  CBX(CN1,CN2,X) 

C»>  EQ.  4 

CBX  =  CN1  *  CF(X)  +  CN2 

RETURN 

END 
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